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Phase-locked loop initialization via curve-fitting 



(57) A method and apparatus for accurately estimat- 
ing the carrier frequency offset and the carrier phase off- 
set of a digitally modulated signal using a signal 
processing algorithm to initialize the state variables of a 
Phase-Locked Loop (PLL) is disclosed. A sequence of 
phase values is estimated from a received sequence of 



symbols and the angular effect due to the modulation 
format is removed from the sequence of phase values. 
A curve-fit algorithm based in one embodiment on the 
RLS algorithm is then applied to a sequence of un- 
wrapped phase values to estimate the carrier frequency 
offset and the carrier phase offset. 
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Description 

BACKGROUND OF THE INVENTION 

s [00011 Coherent demodulation of digitally modulated signals requires a receiver to be synchronized to the carrier 
frequency offset and carrier phase offset of the received signal relative to the transmitted signal. If left uncorrected 
the carrier frequency offset at the receiver may rotate the transmitted signal constellation, which introduces errors each 
time a received symbol rotates past the boundary of a decision region. The carrier phase offset may also introduce a 
fixed rotation to the transmitted constellation, which causes errors due to the misalignment of the decision regions at 

10 the receiver relative to the transmitter. 

[0002] Generally, in the prior art, receivers employ a phase-locked loop (PLL) to acquire and track earner frequency 
offsets and carrier phase offsets. During an initialization period, the PLL locks onto the carrier frequency offset and 
carrier phase offset. Following this period, the PLL tracks these two parameters. 

[00031 For continuous transmission systems, the acquisition period has an insignificant effect on data throughput 
is because the acquisition period is only required once during the entire transmission interval. However, in packet-based 
systems, the acquisition period can have a negative impact on data throughput because each packet requires an 
acquisition phase. . 
[00041 For the foregoing reasons there is a need for improved signal processing methods that may quickly and 
accurately estimate the initial values of the state variables of a phase-locked loop, thereby reducing the acquisition 
20 period. The present invention provides methods and apparatus that meet the aforementioned need. 

SUMMARY OF THE INVENTION 

[0005] The present invention provides methods and apparatus for initializing a phase-locked loop using a signal 

25 processing algorithm. . 

[0006] In one embodiment of the invention, a method for estimating carrier frequency offset and carrier phase offset 
is disclosed. The method comprises the steps of (1) estimating phases of a sequence of digitally modulated symbols; 
(2) removing from each of the estimated phases an angle rotation introduced by a modulation format, wherein the 
phase rotation is computed based on a reference symbol; (3) deriving a set of values from the estimated phases after 

30 removal of said angle rotation, wherein the set of values are a function of the carrier frequency and phase offsets to 
be estimated; and (4) processing the set of values to determine estimates of the carrier frequency and phase offsets. 
In this embodiment, the carrier frequency offset and carrier phase offset are used to initialize a Phase-Locked Loop 

[00071 In another embodiment, apparatus is provided for estimating the carrier phase offset and the carrier frequency 
35 offset The apparatus comprises (1) a phase calculator for estimating phases of a sequence of digitally modulated 
symbols- (2) a remove modulation module for removing an angle rotation introduced by a modulation format to generate 
a sequence of phase values representative of the carrier frequency offset and the carrier phase offset; and (3) an 
estimation module for estimating the carrier frequency offset and the carrier phase offset, whereby the estimation 
module applies a curve-fitting algorithm to the sequence of phase values to generate a linear function dependent on 
40 the carrier frequency offset and the carrier phase offset. In yet another embodiment, the apparatus further comprises 
an unwrap module for modifying the phase estimates generated by the phase calculator module. 
[0008] These and otherfeatures and objects of the invention will be more fully understood from the following detailed 
description of the preferred embodiments, which should be read in light of the accompanying drawings. 

45 BRIEF DESCRIPTION OF THE DRAWINGS 

[0009] 

Figure 1 is a block diagram of a prior art second-order Phase Lock Loop (PLL); 
so Fig Ur e 2 is a block diagram of a PLL initialization module in accordance with the invention; 

Figures 3A and 3B are a block diagram and flowchart, respectively, illustrating a remove modulation module; 

Figures 4A and 4B are a block diagram and flowchart, respectively, illustrating a phase unwrap module; 

Figure 5 is a table of phase values at the output of each module of the PLL initialization module; and 

Figure 6 illustrates the performance of methods and apparatus in accordance with the present invention in graphical 
55 form. 
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DETAILED DESCRIPTION OF THE INVENTION 
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[0010] In describing a preferred embodiment of the invention illustrate in th« * « ■ 

used for the sake of clarity. However the invention is homSC^^ ♦ ! draw,n g*. specific terminology will be 
it is to be understood thaieach specifTc t^ b * 'T* t0 the Speciflc terms so se,ected ' and 

accomplish a similar purpose. e( * u,va,ents whj <* operate in a similar manner to 

ffplLII^ 96nera ' RgUreS 1 ^ 6 «" P^u.ar, the method and apparatus of 

f£ol^ [RLU 1 00 in block diagram form. A Numerica,- 

of the NCO 120 is control by a com"" gna" 7Z^T^tX!rjZ ? T* f?T ^ 0UtPUt 

sS ^'SS^lS^^rS^^ f> computed by the phase detector 130 from the two 
followed by a low-pass-filter Those o^out is JelltTZ? T 11™ im P ,emented . *>' sample, with a mixer 
methods fo" implanting ^ ^e det^r S^^uST ^ ^ ^ inpUt " ° ther 

^whic^^bVr^^^ 

control characterized by a gain (constant* TLTr^lV^ P f,lter 110 Can be realized as a Proportional 
and integrate (PI) control a" I £ c * e " n^S'lS^ ' ^ 11 ° °" 66 rea "' Zed 38 3 P^iona! 
as well. The PLL 100 of Figure 1 can be assodSed SZ^ZSZt ^'T^'" 9 the ' O0p fi,ter 110 ma * be used 
(n) to the phase of the output signal 9 tra " Sfer fUnCtl0n that ,inks th * P»»" °f the input signal x 

riTcSC^ *i£!5^5^^^ haV '" 9 ■""*-»«»•"» (second-order PLL). 
from their initial values K ^ 

carrier phase offset and the carrier frequencv offs^wLS .nit , , PU " the tW ° State variables mav be tne 

steady-state values during the aTquiZr^hasT ' ^ "* eStimated ' 3nd are then brou 9 ht t0 

reference signal is initialized by an initial carrier f™..™™ 5 ! 6 fre< « uenc y of th e input signal. The 

usmg srgnal processing methods as described below it .accordance with ^Figu e 2 P ^ o) ,SCOmpUted 

[001 7] Figure 2 illustrates a method for initializina the PLL 1 nn tiTd. , I r ' . 

sequence of digitally modulated symbols ana genres the oalV I P ^ .'"rtiahzatlon module 200 takes as input a 
zation module 250 to generate the reference siona^i hi * * ' S P3SSed t0 the PLL and NCO initia »- 

module 210 takes as input a SSJSSK SI ^ , ^' n ^ ,c **'* M 
example, from a symbol \xtn^/"^/d^^TLSSr T*** * Ph3Se 6l(n) of the s y mbo1 *<">• For 
generateaphaseestimate wefua ilti^^lg^™?™' tha Ca,CUlate P hasa ™ d "' a 210 wou.d 
an exam P ,e mapping between symbo, *(„) (c^^ 

expressed in degrees, radians or any other scaled version of degrees orrSnT 1 * b ° 

Kao^=rS 

Cat — b V underiying modu.ation 

module 240. The curve-fit algorithm module 240 ^aoomxim^'t?' accomplished by the curve-fit algorithm 

first order po^nomial function refc^ th6 S6qU6nCe ° f ,inear phase estimates °3<") with a 

[0022] In one embodiment, a recursive least-sauare<? (R\ <z\ • . , 

RLS method is thus used to approximate a* m ^^^ eW ^ mn ^ mMM ^ 

[e 3 (0), e 3 (1) e 3 (N-1)]T other Z^^uc^T?, T y Vmg 88 com PO" ent s the linear phase estimates 

method may also be used 9 thS ,east - mean -q"are (LMS) algorithm or the Kalman filtering 
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[0023] The vector y of linear phase estimates is a vector of noisy data that can be represented as a vector u = [u(0), 

u(1) u(N-1)] T plus a noise vector t = [t(0),t(1) t(N-1)] T (i.e. y = u+t). A component u(n) can be represented by the 

linear equation dependent of the pair (9 0 , / 0 ) : u(")=So + n*T s *360*/o where T s is the sampling period (in seconds), 
9 0 is in degrees, and / 0 is the carrier offset frequency in Hertz. 

[0024] The observation vector y can also be modeled with the equation y=H n w+t where H n is an (Nx2) matrix having 

a component h, at row / equal to [1 /*360*T S ] for i=0 N-1 . w is a matrix of parameters to be estimated and is equal 

to [G 0 frf A least-square solution for w is given by the equation w={H N T.H^.H N T y. The estimate w, which is equal 
to [Go , / 0 ] T . may be computed recursively using the RLS or Kalman filtering algorithm. Such algorithms are well known 
to those skilled in the art. Other methods for computing the estimate may be used as well. 

[0025] The estimation of the pair [§„ , / 0 F will be described in accordance with the RLS method. The basic idea 
behind the RLS method is to sequentially update the least squares estimate as new measurements become available. 
The observation equation up to time n may be rewritten as 



y n -i 



• w+u. 



where c n = [1 n*360*Tsp\ It can be shown that that the Least Squares (LS) estimate of w at time n (denoted by w n ) 
can be computed recursively using the following recursions: 



1.y n " 1 = 1 +c n Tp n . 1 c n 

3- P n =Pn-1-YnPo1 C n C/I 

4. OITn-yn-CNTWn., 

5. w n = w rv1 + K n err n 



Assuming P 0 is known or can be accurately estimated and since c n is deterministic, K N can be pre-computed. One 
approach to correctly initialize P 0 is to set it equal to a large diagonal matrix. As an example, letting 



le6 Ol 
0 le6 



and T =1u.s, it can be shown that the first ten values of the 2x1 vector K are given by 



K 0 =[ 0.999000999001 00, 
K 1 =[ 0.00000000770834, 
K 2 =[-0.1 6652789131 688, 
K 3 =[-0.1 998600971 7343, 
K 4 =[-0.1 99880071 64837, 
K 5 =[-0.1 9037646934006, 
=[-0.1 7848855883573, 
K 7 =[-0.1 66597251 1 3632 , 
K 8 =[-0.15549681232540, 
K g =[-0.14540431489250, 



O.OOOOOOOOOOOOOOf 
0.00277777773493] T 
0.001 38865759479] T 
0.0008331 66781 92] T 
0.00055544451 064] T 
0.00039674985045] T 
0.00029756592596] T 
0.000231 44291 728] T 
0.0001 851 5638948f 
0.0001 51 4931 2056] T 



[0026] By applying the vector K to recursions 4 and 5, the estimate w n can be obtained and thus, the pair [e 0 ,/ 0 ] T 
can be computed. 

[0027] The state variables initialization module represented as PLL and NCO initialization module 250 uses the 
carrier frequency and phase offset estimates obtained from the curve fit algorithm module 240 to initialize the state 
variables of the PLL 100. 
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memory elements 301 a, 301 b and Tic The oZZnrtT ^ °° COnfl 9 ured t0 communicate with one or more 
be readily understood from the flowed sn^Tin ^a?SS^^^ 3 °° ^ RgUre 3A °" 
received symbol, the angie rotation relativ tol ^ prevSus SV mhrI J! m ° du,at,on ,0 9 ic 300 estimates, for every 
ments 301a, 301b and 301c are used to ho to vSu« ' ' f „ * J° by thechannel distortion. The memory ele- 

which are updated and ca rriedlZl ^ orevSus ^bni ^ 33 AccumRotation . FirstTimeFlag and State 

current symbol processingTnte^ P pr0Cess,n 9 ,nterva » *V the remove modulation logic 300 to a 

E2J2^^ Binar y Shi« Keying (BPSK) modu.ation 

( M PSK) or Multiple Quadrature Amp't^ 

logic 300 may be applied if the carrier frequency offset / 0 satisNes7<i /Mr ^ l ™ P r ,ementat,on ' tr f remove modulation 
phase shift from one symbol to another will exceed 90 StamJTlSS ^ ^ Camer 1re <* uenc Y the 
Phase shift is due totheLnsmitted sy^ make the determination as to whether the 

effect may be removed if the carrier frequency oXe^ 
,e iP re.ousequation,atrainingsequeL 

time the module is called (step 305) If the n*nZSll 22 Z ■' P erf °™ed to determine if this is the first 

remove modu.ation modu/e 220 se * t^ ou^ phase iu^LlZT °7i * V/f" ^ ' S ^ a " d th * 
the AccumRotation variable and the FirstT^eRa'g va'abTe a^ Let *** ^ ^ * 315 320 

f^^ 



if 0 a fn; > ISO 

elseif $2(11) < -180 

0i(n) = 0 2 (n) + 36-0/ 
else 

9i(n) = ^Cn;,- 

end 



TempPhase = 6 t (n)- State. 

[0033] At step 345, the remove modulation logic 300 ensures that the Dhase different t ^ 
in the range [-1 80» , 1 80°] by using for examole th* a hrZ Z'hk, h £ d,fference - •* - TempPhase, is confined 
at step 350 and the A ccumRoS is updated ^SZEFSZU^. 1 ™ ^ ^ ^ * Ca ' CUlated 
software means, such as the following MATLAB co<S " 355 ™ y be accom P'*hed using a 
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if abs (TempPhase) > 90 

if (0 2 (n) >= 0.0 ) & (02(11) <= 180) 
9 2 (n) « 0i (n) - 180; 

AccumRotation = AccumRotation - 180; 

elseif (0 2 (n) < 0.0) & (0i(n) >= -180) 
0 2 (n) = 0! (n) + 180; 

AccumRotation = AccumRotation + 180; 
end 

else 

02 (n) = 0i (n) ; 

AccumRotation = AccumRotation + 0.0; 

end 



[0034] At step 357, the AccumRotation is confined in the range [-180°, 180°] by similar software means to that 
described above. 

[0035] As illustrated in Figure 3B ( before the remove modulation logic 300 terminates its operation at step 304, the 
40 present state variable is updated at step 360 by setting its value equal to the output phase. 

[0036] Figure 4A illustrates the unwrap phase module 230 of Figure 2. The unwrap phase module 230 contains an 
unwrap phase logic 400 whose operation is described by the flowchart shown in Figure 4B. The unwrap phase logic 
400 performs a modulo operation on the sequence of phase values received from the remove modulation module 220. 
As described by the flowchart of Figure 4B, a bounded phase value 9 2 (n) between [-180°, 180°] is converted to its 
45 absolute value which allows for a first-order polynomial representation of the sequence of phase values. 

[0037] Figure 4B illustrates a flowchart which can be used to implement the unwrap phase logic 400. The test step 
405 determines if the FirstTimeFlag is equal to one, in which case the "Yes" branch is taken. If the FirstTimeFlag is 
equal to zero, the unwrap phase logic 400 proceeds along the "No" branch. 

[0038] If the "yes" branch is taken, at step 41 0 the output phase is set equal to the input phase and at steps 415 and 
so 420, the AccumRotation and the FirstTimeFlag are reset to zero. 

[0039] If the "No" branch is taken, at step 430 a rotated phase value is computed by adding to the input phase value 
the value of AccumRotation. At step 435 a phase difference between the rotated phase value and the present state 
value is computed. Step 430 and step 435 may be accomplished using a software means such as the following M ATLAB 
code: 



55 



B 2 (n) = Q 2 ( n ) + AccumRotation; DeltaTheta = e 2 (n)- State; 
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if f aisrDeltaTheta + 360) < abs (DeltaTheta) ) 
03 (n) = 0 2 (n) + 360; 

AccumRotation = AccumRotation + 350/ 
elseif f abs reel taTheta - < abs (DeltaTheta) ) 

03 (n) = ftfo; - 360/ 

AccumRotation = Accumieotation - 360; 
is else 

03 (n) = 02 (n) ; 

AccumRotation = AccumRotation + 0.0; 

20 end 
end 
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45 



mnlll r ^ PreSSnt State value is updated and the P^ess ends at step 404 

... "constellation points" (column^) ^S^^SSSSSi ( °° " S ? 2) 38 We " as the Emitted symbois. i. 
generated by the phase module 21 0 T^r^ cl^O^ZZ'l J 0 '™ 5 ° 8 C ° nt3 ' nS the phase values 
modulation module 220 of Figure 2. Columbia <xZZ thl JS? . P V3 ' UeS 9 enerated b V th * "move 
230 of Figure 2. " t3,nS the phase values generated by the unwrap phase module 

expresseVlinearrv as a f unctTof ttffCS^SEtf ^r'^ ^ V " U> ^ h 
/o=1/(8T s ) as previously described The noisv ohLe ri?t™L Tf u and tne oarner P hase offset estimate 

represents the'straighMine appnVximatio of £ ^nois^daTa" Tn^ h^, T*?£Z ^ LS 3nd RLS 
obtained by using the RLS method * ' th ® LS and RLS predictor is the Predicted value 

Eg th^^ -ues of the state variables of the PLL a.iows 

offset. The method of the present invent?," rn^be mptemeted n VnuZ^Tn^ T* *"* ^ ph3Se 
hardware or a combination thereof. It can be implemenTed EtJSJlL 2 W3yS inClUdin9 software - 

or implemented as an Appiication Specific Sh 3 Di9it3 ' Signal Proces *or (DSP) 

[0045] Although this invention has been iSSted ^ 

skilled in the art that various changes anTmoSESSS 2 Z Zh TL ,C i embodiments > * wi » *>* «PParent to those 
-"nlsint^^^ 



Claims 



^ 1 . A method for estimating earner frequency and phase offsets of a digita. V modu.ated signai, comprising: 
(a) estimating one or more phases of a sequence of digitally modulated symbols- 

<d> proc ess ,„ g said values ,„ d ., em ,,„ e estimates of , h . ^ « ^ ^ 

2. The method of claim 1 further comprising: 
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(e) initializing the parameters of a Phase-Locked Loop with the estimated carrier frequency and the phase 
offsets. 

3. The method of claim 1 or claim 2 wherein step (c) uses an unwrap phase function to derive said set of- values. 

4. The method of claim 1 or claim 2 or claim 3 wherein the processing of step (d) uses an estimation algorithm based 
on the recursive least-squares method. 

5. The method of claim 1 or claim 2 or claim 3 wherein the processing of step (d) uses an estimation algorithm based 
on the Kalman filtering method. 

6. The method of claim 1 or claim 2 or claim 3 wherein the processing of step (d) uses an estimation algorithm based 
on the least-mean squares method. 

7. A method for recovering a carrier phase offset and a carrier frequency offset comprising the steps of: 

(a) receiving a modulated signal containing a plurality of symbols; 

(b) determining an angular location of first symbol; 

(c) determining an angular location of a second symbol; 

(d) removing the modulation from the second symbol to produce an unmodulated angular sequence; and 

(e) estimating the carrier phase and frequency offsets by curve fitting the unmodulated angular sequence. 

8. The method described in claim 7 further comprising the step of: 

(f) unwrapping the unmodulated angular sequence to compensate for phase discontinuities. 

9. The method of claim 7 wherein step (f) is based on the recursive least-squares method to perform the curve-fitting. 

10. The method of claim 7 wherein step (f) is based on the Kalman filtering method to perform the curve-fitting. 

11. The method of claim 7 wherein step (f) is based on the least-mean squares method to perform the curve-fitting. 

12. The method of claim 7 wherein the carrier phase offset and the carrier frequency offset are used as initialization 
parameters in a phase-locked loop. 

13. An apparatus for estimating a carrier phase offset and a carrier frequency offset, comprising: 

(a) a phase calculator for estimating phases of a sequence of digitally modulated symbols; 

(b) a remove modulation module for removing an angle rotation introduced by a modulation format to generate 
a sequence of phase values representative of the carrier frequency offset and the carrier phase offset; and 

(c) an estimation module for estimating the carrier frequency offset and the carrier phase offset, whereby the 
estimation module applies a curve-fitting algorithm to the sequence of phase values to generate a linear func- 
tion dependent of the carrier frequency offset and the carrier phase offset. 

14. The apparatus of claim 13 further comprising: 

(d) an unwrap module for converting the phase estimates generated by the phase calculator module into 
absolute values. 

15. The apparatus of claim 13 further coupled to a phase-locked loop to initialize the phase-locked loop with the esti- 
mates of the carrier frequency offset and the carrier phase offset. 

16. The apparatus of claim 13 wherein the estimation module applies a curve-fitting algorithm to the sequence of 
phase values to estimate the carrier frequency offset and the carrier phase offset. 
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(57) A method and apparatus for accurately estimat- 
ing the carrier frequency offset and the carrier phase off- 
set of a digitally modulated signal using a signal process- 
ing algorithm to initialize the state variables of a Phase- 
Locked Loop (PLL) is disclosed. A sequence of phase 
values is estimated from a received sequence of symbols 



and the angular effect due to the modulation format is 
removed from the sequence of phase values. A curve-fit 
algorithm based in one embodiment on the RLS algo- 
rithm is then applied to a sequence of unwrapped phase 
values to estimate the carrier frequency offset and the 
carrier phase offset. 
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